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Abstract 

Based on perturbation theory, we study the dynamics of how dark matter and dark energy in the collapsing system approach 
dynamical equilibrium when they are in interaction. We find that the interaction between dark sectors cannot ensure the dark 
energy to fully cluster along with dark matter. When dark energy does not trace dark matter, we present a new treatment 
on studying the structure formation in the spherical collapsing system. Furthermore we examine the cluster number counts 
dependence on the interaction between dark sectors and analyze how dark energy inhomogeneities affect cluster abundances. It 
is shown that cluster number counts can provide specific signature of dark sectors interaction and dark energy inhomogeneities. 

PACS numbers: 98.80.Cq, 98.80.-k 
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I. INTRODUCTION 



There has been convincing evidence indicating that our universe is experiencing an accelerated expansion^!]. In the 
framework of General Relativity, the impart of the acceleration of the universe is given by a kind of mysterious energy 
called dark energy(DE). One leading candidate of such DE is the cosmological constant. However, it is difficult to 
understand such a cosmological constant in terms of fundamental physics. Its value from observation is far below that 
calculated in quantum field theory, what is referred to as the cosmological constant problem. Moreover, using the 
cosmological constant to explain the DE, one is unavoidably led to the coincidence problem, namely, why the vacuum 
and matter energy densities are precisely of the same order today. 

Considering that DE contributes a significant fraction of the content of the universe, it is natural in the framework 
of field theory to consider its interactions with the remaining fields, as e.g., the other large fraction of the universe, 
the dark matter (DM). The possibility that DE and DM interact with each other has been widely discussed recently 
[2]- [23]. It has been shown that an appropriate interaction between DE and DM can provide a mechanism to alleviate 
the coincidence problem [2-6]. However, it has been suspected that the appearance of the coupling between dark 
sectors may lead to the curvature perturbation instability [7] . This problem was later clarified in [8] where it has been 
shown that the stability of curvature perturbation holds while appropriately choosing the forms of the interaction 
between dark sectors and the equations of state (EOS) of DE. Another possibility to cure the instability was suggested 
in [9]. Observational signatures on the dark sector mutual interaction have been found in the probes of the cosmic 
expansion history by using the WMAP, SNIa, BAO and SDSS data etc [11]-[17]. Interestingly it was disclosed that 
the late ISW effect has the unique ability to give insight into the coupling between dark sectors [17]. The small 
positive coupling indicates that there is energy transfer from DE to DM, what can help to alleviate the coincidence 
problem [17, 18]. 

To further unveil the nature of the interaction between DE and DM, we require complementary probes in addition 
to investigating the expansion history of the universe. It has been observed that the growth of cosmic structure can 
be influenced by a nonminimal interaction between dark sectors, what leaves a clear change in the growth index due 
to the coupling [19, 21]. Furthermore, it was suggested that the dynamical equilibrium of collapsed structures such as 
clusters would acquire a modification due to the coupling between DE and DM [22, 23]. Comparing the naive virial 
masses of a large sample of clusters with their masses estimated by X-ray and by weak lensing data, small positive 
coupling has been tightly constrained [23], what agrees with the results given in [17] from CMB. 

The redshift dependence of cluster number counts is another promising tool to discriminate different DE models. 
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Discussions can be found in [24]- [32]. Similar investigation to test the particular scenario of coupled quintessence 
model by using the cluster number counts was investigated in [34-36] , where the interaction between quintessence and 
DM was specifically described by some product of the densities of DE and DM. The sign they used for the coupling 
between dark sectors just describes the energy decay from DM to DE and they concentrated on the inhomogeneities 
in the DM fluid propagating in the scalar field and affecting the scalar field evolution. The cluster number counts 
dependence on the amount of DM coupled to DE was analyzed therein. 

In this paper we will investigate the possibility of employing measurements from cluster number counts to grasp 
the signature of interactions between DE and DM in some commonly discussed phenomenological models. These 
models were proved to have stable curvature perturbations [8] and have been tested by confronting CMB [17] and 
virial mass of cluster observations [23]. In our study we will use the spherical collapse model and consider either DE 
being distributed homogeneously or inhomogeneously. We will discuss the dynamics of DM and the inhomogeneous 
DE to virialization. When DE is homogeneous, DM flows to equilibrium in agreement with that discussed in [38] if 
there is no interaction between dark sectors. If the DE is inhomogeneous, we will show that its dynamics is different 
from that of DM in the virialization in the linear perturbation level. This indicates that DE does not cluster fully 
along with the DM and thus energy is not strictly conserved inside the collapsing system in the linear perturbation 
theory Our result is a progress in the linear perturbation formalism, which has not been reported before. 

In the following discussion of this work, we will use the Press and Schcchtcr formalism [39] to predict the number 
density of collapsed objects. We will develop a new treatment on studying how the structure is formed when DE does 
not trace DM. The Press and Schechter formalism is crude if compared with N-body simulations, while here we do 
not seek the precise confrontations with the observational data, but the evolution of the mass function of collapsed 
objects predicted by the formalism will be helpful for us to understand the influence of the interaction between dark 
sectors and the inhomogeneous DE distribution on cluster number counts. 

This paper is organized as follows: In Section II, we review the formalism of the perturbation theory in the presence 
of interaction between dark sectors. In Section III, we study the dynamics to describe the flow of the collapsing object 
to virialization. In Section IV, we introduce the spherical collapse model in the presence of coupling between DE and 
DM and the Press-Schechter formalism for the galaxy number counts. In Section V, we show the numerical results 
on how the interaction and inhomogeneous DE influence the cluster number counts. In the last section we present 
our summary and discussions. 
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II. PERTURBATION THEORY WHEN DE INTERACTS WITH DM 



In this section, we go over the main results on the linear perturbation theory when DE interacts with DM. The 
detailed descriptions can be found in [8, 19]. 

In the spatially flat Fricdmann-Robcrtson-Walkcr(FRW) background, if there is interaction between DE and DM, 
neither of them can evolve independently. The (non) conservation equation are 



Pr, 



P ' d + m{l + w) Pd = a 2 Q° d , 



(1) 



where the subscript "m" denotes the cold DM, "d" denotes the DE and Q describes the interaction between them. 
Since we know neither the physics of DM nor that of DE at the present moment, we cannot write out the precise form 
of the interaction between them from first principles (see [33] for recent attempts). One has to specify the interaction 
either from the outset [3], or determine it from phenomenological requirements [5, 18]. For the sake of generality, we 
consider the phenomenological description of the interaction between DE and DM in the comoving frame [8, 19] 
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where % is the global expansion rate in the conformal time, £1,62 are small dimensionless constants and T is the 
transpose of the matrix. 

The perturbation equations for DM and DE in the subhorizon scale have been derived in [19], having the form 
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(4) 



l + W 1 + W pd ' " 1 + w ' " p d (1 + w)p d ' 

where the prime denotes the derivative with respect to the conformal time and the gauge- invariant quantities A m , A^ 
represent the density contrast A m ps 5p m /p m = S m , A d w Sp d /p d = S d . The perturbed gauge invariant couplings are 
expressed as [8, 19], 

a 2 SO aI 



Pd 

where r = p m /p d . 

It is useful to rewrite these equations in the real space, 

A' m + V s • V m = 3H6 (A d - A m )/r , 

v^ + nv m = -w^-m^i + ^/r)v m ; (5) 

A^ + (1 + w)V x ■ V d = 3U{w - C 2 e )A d + 3H£ir(A d - A m ) , 

V' d + %V d = -V a v& - -^-V x A d - -^—V d + m{(w- C 2 a ) + 1 ±£p C ° far + 6)1 Vd ; (6) 
1 + w 1 + w { 1 + w J 

where x refers to the conformal coordinates. 

Defining a m = 5p m , a d = 8p d , and assuming that the EOS of DE is constant w' = 0, we can change Eqs.( 5, 6) 

into, 

&m + 3Ha m + V x {p m V m ) = 3Hfaa m + &cr d ), 

j t (aV m ) = -V a (o*) - 3Hfa + &/r)(aV m ); (7) 



& d + 3H(1 + Ci)a d + (1 + w)V x { Pd V d ) = -3Hfaa m + &a d ), 



2 



-(o^) = -V a (otf) - r-^Vx ■ (oAd) + 3ff 
at 1 + w 



2 



(™-e»)+ 1+ 1 ™ u> c ° (ftr + 6) 



(aVi); (8) 



where V x = j;V s . The dot denotes the derivative with respect to the cosmic time, <I> indicates the peculiar potential, 
which can be decomposed into "J = i/} m + ip d , satisfying the Poisson equation [26], 

V^ A = 4ttG(1 + 3w x )a x , (9) 

where a\ represents the inhomogeneous fluctuation field and the subscript "A" denotes DM or DE, respectively. We 
have included the correction from General Relativity. In a homogeneous and isotropic background < ip\ >= 0, since 
< o\ >— 0. For DE and DM, their peculiar potentials read [2G] 



V m = -4^7 / dV- ^— , (10) 

| x — x' | 

fc = -4. G ,^(|±|^i. (11) 



III. DERIVATION OF THE LAYZER-IRVINE EQUATION 



In this section we derive the Layer-Irvine equation [38] when there is an interaction between DE and DM. Layzer- 
Irvinc equation describes how a collapsing system reaches a state of dynamical equilibrium in an expanding universe. 
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In the presence of the interaction between DE and DM, the Layer-Irvine equation can tell us how the coupling between 
DE and DM affects the process of attaining equilibrium as well as the respective configuration. Our derivation will 
be based on the perturbation theory given in Sec. II. 

For DM, the rate of change of the peculiar velocity is given by Eqs.(7) 

^-(aV m ) = -V x {ai> m + a^ d ) - + &M(aV m ). (12) 

ot 

Neglecting the influence of DE and the couplings, the above equation is nothing but the rate of change of the peculiar 
velocity of the DM particle in the expanding universe described by the Newton's law which was the starting point in 
[38]. To derive the energy equation for local irregularities, we follow the same procedure as [38] by forming a product 
of Eqs.(12) with aV m p m i and integrating over the volume. Here i = a 3 dx A dy A dz is the volume element which 
satisfies Jj£ = 3He. Considering Eqs.(l), the LHS of Eq.(12) can be multiplied by aV m and integrated to yield 

-| (a 2 T m ) - a 2 3F(a + fc/r)T m . (13) 

where T m = | J V^pmi is the kinetic energy of DM associated with peculiar motions of DM particles. 

The RHS of Eqs.( 12) can be treated by the same token. Using partial integration, the potential part can be 
changed into 



aV m ^ x (atp m + aip d )p m § = a J V x (p m V m )tp m i + a J V x (p m V m )i() d e. 

Taking account of the first equation in (7), it can turn into 

- J aV m V x (aip m + aip d )p m i = -a 2 (U mm + HU mm ) - a 2 J ip d ^{(r m s) 

+ 3a 2 H{£ 1 U md + &U dm + 2Z 1 U mm + 2&U dd } (14) 

where U mm = \ J cr m ip m i, U dm = J cr d 4' m i , U md = J <r m ip d i, and U dd = \$ <r d ip d e. 
The second term in the RHS of Eq.(12) can be changed into, 

- f(aV m ) 2 3H(^ + &/r)p m e = -a 2 6H(^ + &/r)T m (15) 

Combining Eqs. (13,14,15), we obtain 

f d 

T m + U rnm + H(2T m + U mm ) =- i> dl -(a m i) - 3Hfa + &/r)T m 



dV 

3H {^U md + &U dm + 2^U mm + 2&U dd } . (16) 
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This equation describes how DM reaches dynamical equilibrium in the collapsing system in the expanding universe. 
If the DE is distributed homogeneously, a d = 0, this equation reduces to 

f m + U m m + H(2T m + U mm ) = -3H(6 + 6 A)T m + 6H^U mm . (17) 

For a system in equilibrium, T m = U mm = 0, we get the virial condition. If we take £i = 3£i,£2 = 3^2, we have the 
virial condition (2 + £~i + ^/r)T m + (1 — 2£i)U mm = [23]. Neglecting the interaction £i = £2 = 0, we recover the 
usual virial condition obtained in [38]. One can see that the presence of the coupling between DE and DM changes 
both the time required by the system to reach equilibrium and the equilibrium configuration itself [23]. 

Now we consider the case that DE is no longer homogeneous and has perturbation. The rate of change of the 
peculiar velocity of DE is described by the second equation in (8). Multiplying both sides of this equation by aVd.pdS 
and integrating over the volume, on the LHS we have 



dt 

On the RHS the first term reads, 



-|(a 2 T d ) +3a 2 H{w + t 1 r + e 2 )T d . (18) 



a 2 a 2 f d 

aV d W x (aip m + aip d )pd£ = - — — {U d d + HU d d) - -r- — / ipmsii^ds) 

1 + w 1 + w J at 



■iH{2(C 2 +^)Udd + 2^U mm + ^U md + (C 2 +£ 2 )U dm } . (19) 



1 + w 

For the remaining terms, we have 

- 7 -£- / aV d V x (aA d ) Pd e + 3H \(w - c 2 a ) + 1 + ™~ ° a (ftr + &)1 [{aV d f Pd e 
1 + w J [ l + w J J 

= --^a 2 V d V x (a d )e + 6a 2 H (w-C 2 )+ ^ W ° a for + fr) T d . (20) 
1 + w J l + w 

Combining Eqs.(18,19,20), we arrive at 



(1 + w)f d + Udd + #[2(1 + w)T d + U dd ] = -3H {2{C 2 e + ^)U M + 2^U mm + ^U md + (C 2 + i 2 )U dm ) 

^m-g t (a d s) - c 2 e J V d \7 x (a d )e + 3H[(l + w)(w~2C 2 ) + (l + w~2C 2 )(Z 1 r + &)]T d , (21) 

describing how the DE reaches dynamical equilibrium in the collapsing system in an expanding universe if it partici- 
pates in the structure formation. 

Looking at equations (16) and (21), we see that in non-interacting case (£1 = £2 = 0), when C 2 = 0, DE would 
cluster just like cold DM. Examples of DE models with this property were investigated in [37]. In interacting case, we 
see that for a collapsing system to reach dynamical equilibrium, the time and dynamics required by DE and DM to 



reach equilibrium are different. This shows that in the collapsing system, DE does not fully cluster along with DM. 
Thus, the energy conservation breaks down inside the collapsing system. It would be fair to say that this result was 
obtained in the linear level. It would be more interesting to examine the result in the non-linear perturbation and 
obtain clearer dynamics on how DE participates in the structure formation. 

IV. SPHERICAL COLLAPSE MODEL AND THE PRESS-SCHECHTER FORMALISM 

The simplest analytical tool to study the structure formation is the spherical collapse model. The continuity 
equation for the background DM and DE energy densities are given by 

p m + 3Hp m = 3H(£ip m + ( 2 Pd), 

p d + 31/(1 + w)p d = -3H(Z lPm + &Pd) • (22) 

Considering now the spherically symmetric region of radius R and with the energy density p^ uster = p\ + a\, where 
"A" denotes DM and DE respectively, it will eventually collapse from its gravitational pull provided that a\ > 0. The 
equation of motion for the collapsing model is governed by Raychaudhuri's equation. 

6 = -\e 2 - AnGj2(Px + Spx) (23) 
6 A 

where 9 — 3|| and Raychaudhuri's equation can be rewritten as, 

£ = ^-$>A + 3p A ) (24) 

where R is the local expansion scale factor which is determined by the total energy density in the spherical model. If 
the matter distribution is homogenous in the interior of the spherical region, the local expansion scale factor R will 
have the same value inside the spherical region and from the Birkhoff theorem, the behavior of R is only determined 
by the interior matter and is not affected by the matter distributed outside the spherical region. 

In the spherically symmetric region, supposing the DE distribution is homogeneous (a d = 0) and its evolution is 
still described in (22), the evolution of the DM energy density in the spherical region behaves as 

•cluster I ou ^cluster ou(C ^cluster L £ n \ 

where h = R/R. 

The Raychaudhuri's equation applied to the spherical region has the form 

R = -^G[\pt uster + (\ + w) Pd ]R . (26) 
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where pd is back ground DE energy density. Converting the time derivative into the derivative with respect to the 
scale factor a, we can write R = («) 2 ^r + o4? an d change the Raychaudhuri's equation into 

2a 2 (l + -)R" -R [1 + (3w + l)/r]a = -R[(3w + l)/r + 6 , (27) 
r 

where ( = p^ ster / p m . Therefore, we have 

C' = -(l-6/r)C-3f C+-6A , (28) 
a R a 

where the prime denotes the derivative with respective to the scale factor a of the universe. 

In order to solve these equations, we set the initial conditions R ~ a and R' = 1 at the starting point z = 3200 which 
approximately corresponds to the epoch of matter-radiation equality. At the starting point the cluster is comoving 
with the background expansion. In company with the spherical model, we investigate the linear perturbation in the 
subhorizon approximation satisfying [19] 

d 2 ln5 m 1 3 .,dlnS m ( dlnS,, 



dlna 2 2 2 dlna \ dlna 

- (3^1+6—)— 6 + 3^i6 + 3^ 2 /r-^ 2 — — + l)] + -n m (29) 

r dlna r dlna dlna 2 

and combine with (27,28) by starting from the initial condition Q = 5 m i + 1 in Eqs. (28) till the point when the 
spherical model collapses R(a co ii) ~ to obtain the critical over density S c above which objects collapse. This linear 
extrapolated density threshold is useful in the following computation on the halo abundances. 

If the DE distribution is not homogeneous, the energy content inside the spherical region consists of multi-fluids. 
When DE does not fully trace DM, the four velocity of DE and DM are different ufy ^ «" m ) and 

«(«o = -y(u&,) + «3) ( 3 °) 

where 7 = (1 — v^)^ 1 / 2 is Lorentz-boost factor and v% is the relative velocity of DE fluid observed by the observer 
rested on the DM frame. If DE fully follows DM, v% = 0. Now we consider the non-comoving perfect fluids [20], 

1 (m) — Pm M (m) U (m) 

T$ = p d u? d) u b (d) +p d hf d) (31) 

where 

h ab = g ab + u a u b ig 

the projection operator [20]. Inserting (30) into the second equation of (31) and denoting 
it? n by u a , the energy momentum tensor for DE reads, 

1$ = Pd u a u b + Pd h ab + 2u a q\ d) (32) 



where q a = (p d + Pd)v d is the energy-flux of DE observed by the observer rested in the DM frame. In the above 
equation, we have already neglected the second order terms of v% and assumed that the energy-flux velocity is far 
less than the speed of light v d -C 1,7 ~ 1. In the spherical model, we define the top-hat radius as the radius of 
the boundary of DM. When DE does not trace DM, DE will not be bounded inside the top-hat radius and will get 
out of the spherical region. We assume that the leakage of DE is still spherically symmetric and from the Birkhoff 
theorem, in the spherical region, the Raychaudhuri's equation takes the same form as Eqs. (24). For the energy 
density conservation law, we have, 

VaT$ = Q\x) (33) 

where Q b is the coupling vector as illustrated in Eqs. 2 and "A" denotes DE and DM respectively. The timelike parts 
of the above equation UbQ^x) gi vc 

^cluster + ihp cluster = ^ ^ fluster + £ ap «*«ter) 

p d luster + 3h{l + w)p c d luster = -•&{! + w)p c d luster -2,H{i lP t uster +&p c d luster ) (34) 

where $ — V x v d - The external term incorporating d in the DE density evolution indicates the energy loss caused by 
the leakage of DE out of the spherical region. 

For the spacelike part, only DE has non-zero spatial component and h% Q b d gives 

9fo+4fofo=0 (35) 

where qfy is the dark energy-flux. Assuming that the energy and pressure are distributed homogenously, we obtain, 

& + h(l-3w)d = 3H(Z 1 T + &)d (36) 

where T = p'£ aster / p d uster and we have already used Eqs. (34) and kept linear order terms of From Eqs. (36), if i9 
vanishes initially, 1? will keep to be zero at all the time during the evolution, DE will fully trace DM. However, in most 
cases, even in linear region there is a small difference between v d and v m that the initial condition for ■& is non-zero. 
We take the initial conditions for # as $ ~ k(vd — v m ) ~ — <5 rra i/200 < 0, which is obtained from the prediction of 
linear equation with k = IMpc" 1 at z% — 3200. ■& is much smaller than 1, \&\ << 1, and the negative value of d means 
that at the initial moment DM expanded faster than that of DE. 

Defining ( m = p^ ster / p m , Cd = p d uster / ' p d and converting the time derivative from ^ to ^ we have the evolution 
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of DM and DE in the spherical region described by 



+ ^(1 -3t«)0 =-(£^ + 6)0, 
K a 

C = - {1 - 6/r} C™ - 3^C m + -6Cd/r, 
a K a 

C' d = -(l + w + Cir)Cd - 3(1 + ^^Crf - -eiCmr - 0(1 + «;)&• 
a K a 



(37) 



When DE fully traces along DM, i? = 0, only the last two equations above are needed to discuss the structure 
formation. The Raychaudhuri's equation now becomes 



2a 2 (l + -)R" -R'[l + (3w + l)/r] a = -R [(3w + l)Q/r + Cm] • 



(38) 



Taking the same initial conditions as above for the spherical region of radius R and in addition adopting the adiabatic 
initial conditions, Sdi = (1 + w)8 m i, and Cmi = b~mi + 1, Qi = Sdi + 1 which lead to Qi = (1 + w)Cmi — w, we can study 
the spherical collapse of DM and DE . However, it is important to note that due to the coupling, Q may get negative 
at sometime during the collapse. To avoid this unphysical point, we make a cut-off of the coupling when 6 becomes 
negative to guarantee Cd >= in our analysis. For linear perturbation, we adopt the equations in the subhorizon 
approximation through [19] 



d 2 lnS 7 



dlnS„ 



dlna 2 \ dlna 

3[exp(ln-^-) - 1] 



1 3 . 



6+366+3£ 2 2 /r + 6( 



dlnS„, , £ 2 > dlnSrr, f? dlnSd , 5h , 

— — - (36 + 6^)—^ + 3^ — J-expQn^) 
dlna r dlna r dlna o m 



dlnH 
dlna 



i)-6 



dh 



dlna 



Q m + (1 - Q m )exp(ln— ) 

Om. 



, (39) 



d 2 ln6d 
dlna 2 



dlnS a 
dlna 



-w(l - fi m ) 



dlnSd , , 1 , .3 

— + (l + w)- 

dlna 2 



36 + 6£ir + 6uj - 3C a + 3(C e - C a — — + 7— ^ — 

1 + u> 1 + w dlna 



£l m exp{ln-P-) + (1 - fi m ) 

dlnSd 
dlna 



3r^iexp(ln 



§m \ dlnS m 
Sd dlna 



3( 



dlna 



dlnH 



dr 



+ l)(w - C e 2 ) + 36(^77- + r + — ) - 3£ 



(-» — + 1 > 
dlna 



k 2 C 2 
a 2 H 2 



dr 



dh 



exptln—) 

Od 



iv - C 2 + Zir(l - exp{ln^)) 

Od 



(1 - 3u,) - 3^-^(1 + u, + f ir + 6) - 3(eir + 6)~ TX 

1 + w 1 + ama 



(40) 



We see that the DE and DM perturbations are entangled. In the subhorizon approximation, it was observed that the 
influence of the DE perturbation is small if compared with that DM [19]. Using ( 39, 40) together with ( 37, 38), we 
can obtain in the spherical model with inhomogeneous DE distribution the linearly extrapolated density threshold, 
S c (z) = 5 m (z = Zcoii), above which structure collapse. 

The spherical collapse model was used by Press and Schechter[39] in providing a formalism to predict the number 
density of collapsed objects. Although this formalism is crude as found in [40] etc., it predicts the evolution of the 
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mass function of collapsed objects well enough for the purpose in this paper to see how the interaction between DE 
and DM influences cluster number counts. 

The comoving number density of collapsed dark halos of mass M in the mass interval dM at a given redshift of 
collapse is given by [39] 



dn(M,z) _ [2 Pm S c ^ S 2 /2a 2 



R da 

7Zr 



(41) 



dM V 7T 3M 2 a 

where ~p m is the comoving mean matter density at a given redshift. In most cases it is a constant and equals to the 
present mean matter density, but this is not true when DE interacts with DM [35]. The quantity a = cr(R, z) here is 
the variance smoothed over radius R. It has an explicit form [41], 

/ R V 7(«) 

^ = H^M^ D[ZI (42) 
where trg is the variance over a sphere with radius R = 8/i _1 Mpc and D(z) is the growth function defined by 

D(z) = 5 m (z)/6 m (0). The index 7 is a function of the mass scale and the shape parameter T of the matter power 

spectrum [41] 

R 



j(R) = (0.3r + 0.2) 



2.92 + log 
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(43) 



^8/i- 1 Mpc / 

We will use T = 0.3 throughout our analysis. The radius R at given M can be calculated by the relation [41], 

1 ( Mh \ 1/3 

where is critical density at present and M s is the solar mass. 

The Press-Schechter formalism presents us the comoving number density of halos, which can be compared with 
astronomical data. In order to do the comparison, we calculate the all sky number of halos per unit of redshift in the 
mass bin 

^ = / dn-^r [ n(M)dM, (45) 



dz J dzdQ _ 

where the comoving volume element per unit redshift is dV/dzdil = r 2 (z)/H(z) and r(z) is the comoving distance 

In the next section, we will present numerical results to see the effect of the interaction between DE and DM. 
In each situation we will study the cases with homogeneous and inhomogeneous DE distributions. When the DE is 
distributed inhomogeneously, it will participate the collapse and the structure formation. We are not going to seek 
precise confrontations with observational data in this work, but to understand the influence of coupling between dark 
sectors on cluster number counts. 
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V. CLUSTER NUMBER COUNTS 



In this section we present numerical results on the cluster number counts due to different interactions between 
DE and DM. We limit our attention to three commonly used forms of interaction between dark sectors, with the 
interaction proportional to the energy densities of DE, DM and total dark sectors respectively. In our numerical 
calculation, we consider that all models have the same <jg = 0.8 . In the end of this section we also discuss the results 
by normalizing the halo abundance to the present value. 

A. Interaction proportional to the energy density of DE (£1 =0,^2 / 0) 

When the interaction is proportional to the DE density, it was found that the curvature perturbation is always 
stable for both quintessence and phantom DE EOS [8]. When we consider DE EoS w > — 1, we have the numerical 
results as shown in Fig 1. 

In the case when DE distribution is homogeneous, the results are shown in solid lines. With positive coupling 
(£2 > DE decays to DM), we see in Fig la that the critical mass density decreases compared with the ACDM model. 
This means that with the positive coupling, cluster can be formed easier than that of ACDM model. However when 
the coupling is negative, with DM decays to DE, 5 C is higher than that of ACDM, which means that it is harder for 
the cluster to be formed. The impact of the coupling between dark sectors is significant in causing the difference in 
the cluster number counts. For positive coupling (£2 > DE decays to DM), the cluster number counts are bigger 
than that of ACDM model, which is clearly shown in Fig lb-Fig Id. However when the coupling is negative^ < 
with DM decays to DE), the situation is opposite and consequently cluster number counts are smaller than that of 
ACDM model. 

When the DE is distributed inhomogeneously, we need to consider the DE perturbation and its participation in 
the structure formation. Results are shown in dotted lines in Fig 1. We see that 5 C is slightly larger than the results 
without DE perturbation. The difference brought by the DE inhomogencity is really small if compared with that of 
the interaction. 

For the DE EoS w < —1, the numerical results are shown in Fig 2. When the DE distribution is homogeneous, 
the influence brought by the interaction between dark sectors on the evolution of critical overdensity S c and the 
galaxy number counts agrees well with the situation when w > — 1. However when we consider the DE is distributed 
inhomogeneously, the model with clustering DE with w < — 1 has a bit more number counts than the homogeneous 
DE distribution appeared in the large mass bin, although the difference is very small especially for the case of small 
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Figure 1: The numerical results for the coupling proportional to the energy density of DE when the EOS of DE w > — 1 . 
The solid lines show the results calculated with homogeneous DE distribution and the dotted lines are for the inhomogeneous 
distribution of DE. 



mass bin. 

Comparing with the interaction, we see that the influence brought by the difference between the homogeneous and 
inhomogeneous DE distributions on the critical density S c and the number of galaxy clusters is very small. This in 
fact tells us that the fluctuation of DE field ad is small. In eq.(21), this small ad leads to U m d ~ Vdm ~ Udd ~ Td — > 0, 
which means that the DE plays very little role in the virialization of the structure. The structure formed is almost the 
same as that of the cluster with homogeneous DE distribution. The little effect of the inhomogeneous DE distribution 
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Figure 2: The numerical results for the coupling proportional to the energy density of DE when the EOS of DE w < — 1 . The 
solid lines show the results calculated when DE is distributed homogeneously and the dotted lines are for inhomogeneous DE 
model. 

in the structure formation is consistent with the finding that the DE perturbation is small compared with the DM 
perturbation in the subhorizon approximation [19]. 

B. The interaction proportional to the energy density of DM (£i > 0, £2 = 0) 



In this case, the stable curvature perturbation can only hold when DE EOS is smaller than -1[8]. The coupling 
between DE and DM should be positive to avoid the negative energy density of DE at the early time of the universe[18] 
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The numerical results by choosing this coupling arc shown in Fig 3. 

The solid lines are for the homogeneous distribution of DE. For small coupling, the lines are close to that of ACDM 
result. We observe that more positive coupling leads to more number counts of clusters Fig 3c. This can be understood 
from 8 C and also the ratio 5 c {z)/ a%D(z), which are smaller at high redshift for more positive £1. This result is also 
consistent with the case when the interaction is proportional to the DE energy density. 

The dotted lines are for the situations when DE is distributed inhomogeneously and we need to consider the DE 
perturbation and its participation in the structure formation. From the behavior of S c we see that S c is suppressed at 
low redshift. This suppression is more obvious than the case of homogeneous DE. In the ratio 8 c (z) / asD(z) shown in 
Fig 3b, the ratio due to the inhomogeneous DE is smaller than that of the homogeneous DE. This leads to a larger 
cluster abundance due to the inhomogeneous DE distribution than that of homogeneous DE distribution as shown in 
Fig 3c. 

C. Interaction proportional to the total energy density of dark sectors (£ = £i = £2 > 0) 

Here, again, the curvature perturbation can only be stable when DE EOS is smaller than -1 [8] and the coupling has 
to be positive [18]. The numerical results are shown in Fig 4. The results are similar to the case when the coupling is 
proportional to the energy density of DM. With more positive coupling, the critical overdensity is smaller which leads 
to the bigger galaxy number counts than that of the ACDM model. For the same coupling, if the DE is clustering, the 
critical overdensity S c and S c (z)/a s D(z) will be suppressed at small redshift. Therefore it is observed in Fig 4c, Fig 4d 
that at small the inhomogeneous DE distribution leads to a larger cluster number abundance as compared with the 
homogeneous DE distribution. 

D. Normalization 

So far in the above discussions we have fixed as = 0.8 as given by present observations [42]. Considering that at 
redshift zero all models have the same number density of halos, the normalization is required to be done by adjusting 
us in each model such that 8 c /a% is equal to the fiducial a$ — 0.8 ACDM model 

_ 5c,model{z — 0) _p ,. ■. 

C ,A( Z = 0) 

We dedicate the rest of this section to discuss the implications on our results by normalizing models to the same 
halo abundance at z = 0. It is interestingly observed that the qualitative structure of curves changes dramatically 
depending on the choices of the normalizations, which was also found in studying the number counts in homogeneous 
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Figure 3: The numerical results for the coupling proportional to the energy density of DM. The solid lines shows the results 
calculated with homogeneous DE distribution and the dotted lines are for inhomogeneous DE. 

and inhomogeneous DE models without coupling to DM [27]. 

In our models, we find in Fig 5a that when the coupling is proportional to the energy density of DE(£i = 0, £2 7^ 0), in 
the normalization of local halos abundance, the departure from the cosmological constant model due to the interaction 
between dark sectors becomes smaller especially in the low mass bin when compared with the results by fixing erg- 
The major impact due to the coupling emerges in the massive galaxy cluster as shown in Fig 5b. The signature of 
the DE perturbation is still not obvious which is similar to the result by fixing a$. 

When the coupling is proportional to the energy density of DM, we see from Fig 6 that the departure from the 
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Figure 4: The numerical results for the coupling proportional to the energy density of total dark sectors. The solid lines shows 
the results calculated with homogeneous DE distribution and the dotted lines are for inhomogeneous DE. 



cosmological constant model due to the interaction between DE and DM (solid lines) is suppressed if compared with 
the result by fixing as Fig 3. This shows that fixing the halo abundance, the coupling does not show up in the galaxy 
number counts. However, in this case, the difference caused by the clustering DE (dotted lines) is enhanced with large 
deviations from the cosmological constant model. When the coupling is proportional to the energy density of total 
dark sectors, the qualitative results are the same. 
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Figure 5: The results for the coupling proportional to the energy density of DE by fixing the local abundance. 
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Figure 6: The results for the coupling proportional to the energy density of DM by fixing the local abundance. 
VI. CONCLUSIONS AND DISCUSSIONS 

Assuming first that DE component is homogeneous, we have investigated how DM approaches the dynamical 
equilibrium in the collapsing system. In the presence of coupling between DE and DM, the flow of mass and energy 
between the components changes the time required by the system to reach equilibrium and also the equilibrium 
configuration itself. Our discussion is based on the perturbation theory developed in [8, 19]. We reproduce the 
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Layzer-Irvine equation [38] when the interaction is neglected. In the homogeneous case, the DE component flows 
progressive out of the overdensity which leads the energy conservation law invalid in the collapsing system [43] . It is 
expected in the models of DE interacting with DM, if the DE is distributed inhomogeneously, it may cluster along 
with the DM and thus we can avoid the energy non-conservation problem. To see whether this speculation is correct 
or not, we have further explored the dynamics of the DE that presents inhomogeneitics at cluster scales and found 
that in the presence of the coupling between DE and DM, the time and dynamics required by the DE in the system 
to reach equilibrium is different from that of DM. This shows that in the collapsing system, DE does not fully cluster 
along with DM, even in the presence of the interaction between them. The energy is not strictly conserved in the 
collapsing system. This result is based on the linear perturbation formalism, where the DE perturbation can be 
neglected in the late universe [19]. 

We also developed a novel treatment of the spherical collapse model which consists of multi-fluids with different 
collapsing velocities. When DE docs not trace DM in the spherical collapse, we studied the dependence of the cluster 
number counts on the coupling between DE and DM and on the DE inhomogeneities by using the Press-Schechter 
formalism. In most of our work we have assumed that the models are normalized to have practically the same erg = 0.8 
as given by present day observations. We have shown that when the coupling between dark sectors is proportional to 
the energy density of DE, there is a significant dependence of cluster number counts on the interaction. If the energy 
decays from DE to DM, the number of cluster counts increases. The dependence of cluster number counts on the DE 
inhomogeneities is negligible compared with the interaction influence. This in fact shows that the fluctuation of the 
DE field is small and the inhomogeneous DE plays little role in the virialization of the structure. For the coupling 
between DE and DM proportional to the energy density of DM and total dark sectors, we have observed that more 
energy decay from DE to DM leads more abundance of the cluster number counts compared with the ACDM model. 
In addition, the signature of the DE inhomogeneities is enhanced. This enhancement can be understood from the 
dynamics in which it shows that the DM gravitational effect pulls the inhomogeneous DE to move and collapse along 
with DM. 

We have further studied these models enforcing the normalization to reproduce the same abundance of galaxy 
clusters at redshift zero. The qualitative structure changes in choosing different normalizations [27]. Comparing with 
the interaction effects in fixing as,, we have verified that the departures from the ACDM model are suppressed, the 
coupling between dark sectors while fixing the local abundance usually does not show up in the galaxy number counts, 
except when the interaction is proportional to the energy density of DE in the super mass bin. In fixing the local 
abundance, the influence caused by the inhomogeneous DE is enhanced with large deviations from the cosmological 
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constant model when the coupling is proportional to the energy density of DM and total dark sectors. 

The Press-Schechter formalism we employed is good enough for the purpose of this paper to illustrate how the 
interaction between dark sectors and the DE inhomogeneities influence cluster number counts, but it is still not easy 
to seek precise confrontations with observational data at the present moment. Recently in the uncoupled model, 
the statistical analysis was performed on the sensitivity of cluster counts to DE perturbations by employing a set of 
cosmological parameters using the Fisher matrix method. It was argued that the impact of DE fluctuations is large 
enough, which could already be detected by existing instruments [32]. It is of great interest to extend their analysis 
to examine the interaction between DE and DM effect and the DE inhomogeneities influence. Such a study is under 
way. 
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